Rydberg States of Lithium and Other Systems 
The interesting calculation of retardation corrections to excited atomic state energies is facilitated by 
the existence of good calculations of the nonrelativistic and nonretarded energies, to which 
retardation corrections can be added. Here I describe a perturbation method, applicable when the 
angular momentum is high enough, which generates an asymptotic series for the energy. It makes use 
of generalized polarizabilities of the atomic or ionic core (numerically obtained) but otherwise is 
completely analytic. It is applied here to the lithium atom in detail, and its application to hydrogen 
molecular ions is outlined. 
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I. INTRODUCTION 

An important discussion component of this Topical Group meeting is the interaction of a highly 
excited single atomic electron with the remaining core electrons when retardation effects are 
important. Such a Rydberg state can be investigated by measuring the deviations of fine-structure 
levels from their unperturbed, hydrogenic values; surprisingly high precision has already been 
achieved in helium [1] and lithium [2] atoms. To carry out the theoretical analysis of these states it 
is still necessary to obtain a good solution of the nonretarded problem, to which the retardation 
corrections can be added. Drake [3] has perfected a variational method for finding accurate energies 
of such states, and his results are definitive, at least for the lower-lying levels. I will give a brief 
summaiy here of another method that has been used; it has the advantages of being largely analytic 
rather than numerical and of improving as the Rydberg electron moves further from the core. 

Consider the lithium atom, with two electrons in a Is 2 'S 0 core state and the third electron in 
a highly excited state with quantum numbers ( N,L ), described [4] by the following nonrelativistic 



Hamiltonian (Z=3 for lithium): 
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The first bracket is the Hamiltonian of the core, H„(12), the second of the Rydberg electron, Ho(x) 
and the braces contain the interaction between the two parts, H^ , which decreases rapidly with 
increasing distance of the third electron, x. To a pretty good approximation, the energy of this system 
is e 0 -1/A/ 2 , where e„is the energy of the core, and the third electron moves in an almost pure 
Coulomb potential, 2/3 shielded by the core. Our problem is to find the small, L-dependent 
corrections to this zero-order energy. 


II. THE METHOD 


We first expand the interaction in a multipole series, with x»r 12 : 
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and then calculate perturbations in the usual way, beginning with second order: 
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where v n (x) = <0\H int \n> . In principle, we could evaluate this sum, including excited states of both 
the core ( n ) and the outer electron (N',L ^ . Not only would this be somewhat difficult but it would 
have to be repeated for every value of N and L in which we were interested. Instead, we make the 


assumption (justified after the fact) that the excitation energies of the core are larger than those of 
the outer electronic system. (This looks reasonable, but there are certainly some huge excitation 
energies of the Rydberg electron up into the continuum that we hope are “negligible”). Then we 
expand the denominator of Eq.(3) in a binomial series: 
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In the first term above it is easy to apply closure over the intermediate states N’L ’ to obtain an 
effective “adiabatic” potential 
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whose expectation value (NL\V ^x^L^ gives the leading correction to the energy. The second sum 
in Eq. (5) cannot be carried out to °° because the higher terms would diverge; the series is probably 
asymptotic and should be terminated at its smallest term. 

The second term in Eq. (4) seems harder to convert into effective potential form because of 
the energy-dependent factor, but a simple trick solves this problem: 
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By converting the energy dependence to a commutator we are able to carry out the closure sum; I 
think the first time this trick was used may have been by Mittleman and Watson [5], After 
commutation and integration by parts, this “first nonadiabatic term” takes the following form: 
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where P, is just like a, but with one more power of the energy difference in the denominator. And of 
course higher terms can be computed in a similar way. (This way of describing the process of 
constructing the effective potential is easy to understand, but I prefer a more formal method using 
projection operators that is described in detail in my chapter of the Casimir book [6].) 

In any case, the effective potential can be written out to the term of order x' 8 as follows: 
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In a series of papers [7] which are summarized in perhaps excessive detail in Ref. [6] I applied this 
method to the Rydberg states of helium, and I even included terms up to order x' 10 with considerable 
success. Of course, in that case all the coefficients that appear in the effective potential can be 
evaluated exactly and analytically, since the core is simply a hydrogenic ion. In the case of lithium, 
the core is a two-electron system, and the sums giving the various coefficients must be performed 
numerically. To calculate them, we used the method of pseudostate summation; that is, we replaced 
the complete set of states \ri> by a finite and discrete set of pseudostates. These are in the form of 
prediagonalized Hylleraas expansions, and e n is replaced by the corresponding expectation value of 
the Hamiltonian. This method converges well, and the effective potential is (for 7 Li, with mass- 
dependent effects included) the following [8]: 
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It is now only necessary to evaluate the expectation values of the inverse powers of x in the 
appropriate hydrogenic state (N,L) to obtain our values of the energy shifts, and there are very well 
known formulas for that. (Actually, small second-order terms, given analytically, must also be 
included to the accuracy needed here.) 

After this nonrelativistic calculation it is possible to add some other small corrections [9], 
These include relativistic effects, which we treat with the Breit-Pauli operator, the Lamb shift, and 
finally, the retardation correction, which is the main topic of this meeting. Remarkably, inclusion of 
these three corrections makes a significant improvement in the agreement of theory with experiment 
in the one interval (1 OH- 101) for which the estimated errors are small enough. The results are given 
in Table I. You can see that experiment agrees well with the theory, including the small corrections, 
for the interval 1 OH- 101, while the results are not definitive for the 10G-10H interval. 

m. MOLECULAR HYDROGEN ION 

Inspired by recent experiments on the fine structure of Rydberg states of the hydrogen molecule by 
Lundeen’s group [10] and their analysis in terms of an effective potential similar to that described 
above, we [1 1] have tried to calculate the dipole polarizability of the hydrogen molecular ion, H 2 + . 
We had tried to use our Hylleraas variational method to calculate the ground state of this ion some 
years ago without quantitative success; this type of trial function is clearly designed for atomic 
systems rather than the more localized motion of the very massive nuclei. But we have recently been 
more successful by allowing very high powers of the intemuclear coordinate to be included. For 
example, the function R N e aR , with a=N/2, produces a function peaked at R=2, close to the 
adiabatic equilibrium point, and if N is large enough the function can represent the Gaussian-like 
vibrational wave function of the ion quite well. Using a ground-state wave function with 308 terms 



and including up to 364 intermediate p-pseudostates, we find the polarizability of the ion to be 
3.16795 a 0 3 , to be compared with the experimental result [10] of 3.1681(7) q? . This agreement is 
good, although our result is not yet completely converged. I am not sure whether it is practical to 
hope that in the future theory and experiment will be able to observe Casimir effects in H 2 . 
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TABLE I 


Comparison of level differences for 7 Li, in MHz, between theory and experiment [2], 
The standard deviations of the differences are experimental and theoretical combined 
in quadrature. 


Interval 

Theory 

Experiment 

Experiment-Theory 

10G-10H 

339.697±.115 

339.7186±.0031 

0.02±0.1 1 


1 OH- 101 


109.2137±.001 1 


109.2140±.0047 


0.0003±0.0048 















